Pairing strengths for a two orbital model of the Fe-pnictides 
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Using an RPA approximation, we have calculated the strengths of the singlet and triplet pairing 
interactions which arise from the exchange of spin and orbital fluctuations for a 2-orbital model of 
the Fe-pnictide superconductors. When the system is doped with F, the electron pockets become 
dominant and we find that the strongest pairing occurs in the singlet d-wave pairing and the triplet 
p-wave pairing channels, which compete closely. The pairing structure in the singlet d-wave chan- 
nel corresponds to a superposition of near neighbor intra-orbital singlets with a minus sign phase 
difference between the dxz and dy^ pairs. The leading pairing configuration in the triplet channel 
also involves a nearest neighbor intra-orbital pairing. We find that the strengths of both the singlet 
and triplet pairing grow, with the singlet pairing growing faster, as the onsite Coulomb interaction 
approaches the value where the 5* = 1 particle-hole susceptibility diverges. 

PACS numbers: 71.10.Fd, 71.18.+y, 71.20.-b, 74.20.-z, 74.20.Mn, 74.25. Ha, 75.30.Fv 



Recently, a new class of superconductors involving a 
family of Fe-based oxypnictides has been discoveredlU 
1 i i, i, S 0, i i- With as high as 55is:0, 
the mechanism of superconductivity is Ukely to be elec- 
tronic in origin and consequently, these materials have 
generated tremendous excitement. Moreover, experimen- 
tal results including specific heatflT 
spectroscopy [l^ . high- field resistivity. 



11 1 , point-contact 
li, Q and NMR 



measurements suggest the existence of unconven- 
tional superconductivity in these materials. Further- 
more, transport [lij and neutron scattering 17 1 measure- 
ments in LaOFeAs have shown the evidence of spin- 
density-wave (SDW) magnetic order below T = IZIK . 
An experimental determination of the orbital and spin 
state of the Cooper pairs, however, has not yet been 
made. 

Band structure calculations show that the Fermi sur- 
face of F doped LaOFeAs consists of two nearly concen- 
tric hole cylinders surrounding the F point and two ellip- 
tically distorted electron cylinders around the M point of 
the 2Fe/cell Brillouin zone. Electronic transitions involv- 
ing states on one or between two of these Fermi surface 
sheets lead to q-dependent structure in the spin and or- 
bital susceptibilities. For small doping, the electron and 
the hole fermi surfaces are of comparable sizes, and their 
nesting can give rise to the observed SDW order in the 
undoped material H, 17 1 . Upon further doping, the two 
hole pockets shrink, the electron fermi surfaces become 
dominant llSll . and the system exhibits superconductiv- 
ity. Ref. [l9| suggests that a triplet p-wave pairing state 
is obtained on the electron Fermi surfaces due to the 
ferromagnetic spin fluctuations. Other related possibili- 
ties have also been discussed in the literature, including 
inter-orbital on-site triplet pairing (20l . 21 1 , and a s-wave 



pairing state which changes sign from the electron to the 
hole pockets [22l |. 

Recently, we have introduced a tight-binding model 



[23| with "(ia;z" and " dy^" orbitals on a two-dimensional 
square lattice of "Fe" sites. This simple tight-binding 
model correctly reproduces the topology of both the elec- 
tron and the hole fermi surfaces. It also reproduces the 
van Hove singularities obtained in bandstructure calcula- 
tions. For low doping, when the electron and hole pockets 
are comparable, RPA calculations show enhanced SDW 
fluctuations at the wave vectors (tt, 0) and (0, tt), defined 
in the convention of one Fe atom per unit cell [l^]. In 
this work we investigate the nature of the pairing state 
when this model is further doped. With on-site inter- 
orbital and intra-orbital Coulomb interaction terms, we 
use the RPA approximation to study the effective pair 
interaction vertex induced by the spin and orbital fluc- 
tuations. We find that when the doping is increased and 
the electron pockets become larger, the leading pairing 
instability occurs in the singlet d-wave and the triplet 
p=wave channels. The pairing strength for both chan- 
nels increases as the system approaches an instability in 
the S ~ 1 particle-hole channel, with the singlet d-wave 
channel growing faster than the triplet p-wave channel. 

Model Hamiltonian - Our tight-binding model Hamil- 
tonian describes a square two-dimensional "Fe" lattice 
with two orbitals per site 

k(T 

(1) 

Here a is the spin index, Ti are Pauli matrices and — 
[c?J,^(A:), is a two-component field, which describes 

the two degenerate "(ia;z" and "dyz" orbitals. The ma- 
trix elements of Hq, e+(k) ~ —{ti+t2){cos kx + cosky) — 
At 3 cos kj; cos ky , e_(k) = — (ti — t2){coskx — cosky) and 
£xyO^) — — 4t4 sin kx sin ky are parametrized by four hop- 
ping paramters i^, i = 1, • • • , 4. 

This free fermion Hamiltonian is diagonalized by intro- 
ducing a canonical transformation to the band operators 
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with 



(2) 



6_(fe) 



i" - 



(3) 



the wave-function of the v band with v = ±1, and 71^0-, k 
the annihilation operator of an electron with spin a and 
wave-vector k in the v band. With the inclusion of a 
chemical potential n, the band part of the Hamiltonian 
becomes 



= ± (4) 



with E±{k) = e+{k) ± ^e'i(k) + ely{k). The tight 

binding parameters ti can be adjusted to fit the Fermi 
surface obtained from LDA band structure calculations 
[lil . 22, 24 1 . In this work, we will take the parameters 
ti = —1, t2 = 1.3, — — —0.85 and measure energy in 
units of With a chemical potential n — 1.45, one has 
a filling of 2 electrons per site, a Fermi surface similar to 
bandstructure calculation of lightly doped LaOFeAs and 
a peak in the bare spin susceptibility for q — (tt, 0) and 
(0,7r). Here we will take fi = 2.0 which corresponds to 
having 2.32 electron persite and gives the Fermi surface 
shown in Fig. [T] (a). There are four Fermi pockets in the 
Brillouin zone: ai around (0,0) and a2 around (tt,tt) 
are hole pockets associated with E-{kf) — 0, while (3i 
around (tt, 0) and P2 around (0, tt) are electron pockets 
given by E+[kf) — 0. For this minimal model, we will 
include only onsite intra and inter orbital Coulomb in- 
teractions, which will both be set equal to U and we will 
neglect the Hunds rule coupling. In this case, up to a shift 
of the chemical potential, the interaction can be written 



Hint — 



E 



(5) 



An important feature of this two band model is the 
nontrivial C4 rotation symmetry of the two orbitals. Un- 
der a 90° degree rotation, the two orbitals transform as 
\xz) — > \yz) and \yz) —\xz). Correspondingly, in the 
Hamiltonian, e+(k) has s-wave symmetry and e-(k) and 
exyi}i.) have d-wave symmetry, which together perserve 
the point group symmetry of the Hamiltonian. Con- 
sequently, the wave functions j^. of the energy eigen- 
states also have nontrivial structure in the Brillouin zone. 



which can be determined by the direction of the vector 
n(k) = (e_(k), ea;y(k)). If we consider the orbital degree 
of freedom as a pseudo-spin, the electrons in the lower 
band always have a "pseudo-spin" anti-parallel to n(k). 
For the parameters we are using, the distribution of the 
unit vector n(k) ~ n(k)/|n(k)| is shown in Fig. [T](b). 
For example, at the wavevector k = (7r/2,0) we have 
e_(k) < and e2:y(k) = 0, which means the upper band 
is formed from xz orbitals and the lower band from yz 
orbitals. From Fig. [1] (b) we can see that the electron 
pocket Pi {P2) is formed mainly from xz (yz) orbitals, 
while the hole pockets are formed from "d-wave" super- 
position of the two orbitals. This point will be important 
for understanding the pairing symmetry. Since this non- 
trivial structure of the wave function originates from the 
symmetry of the two orbitals dxz and dyz, we expect it 
to be qualitatively correct even beyond the present two 
orbital model. 
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FIG. 1: (a) The Fermi surface of the 2-orbital model on the 
large IFe/cell BZ. Here, the ai,2 surfaces are hole Fermi pock- 
ets given by E-(kj) — and the /3i,2 surfaces are electron 
Fermi pockets given by Ej^(kf) = 0. In this paper we have 
set ti = -1,£2 = 1.3, is = U = -0.85 and n = 2. (b) Wave 
function distribution in the Brillouin zone. The arrows show 
the direction of the vector (e_(k), exyi}^))- When an arrow is 
pointing up (down) at some k point, the eigenstate of upper 
band iJ+(k) consists of pure xz {yz) orbitals. The Brillouin 
zone is shifted by (7r/2,7r/2) for convenience. 



In the following we first discuss the particle-hole sus- 
ceptibility and calculate it using an RPA approximation. 
Then using the pairing interaction associated with the ex- 
change of these particle-hole excitations, we examine the 
strength of the pairing in the singlet and triplet channels. 

One loop and RPA susceptibilities - Because of the 
two-orbitals, the generic form of the susceptibility de- 
pends on four orbital indices p,q, s,t equal to 1 or 2 for 
dxz and dyz, as well as spin indices: 



dre 



•^.a,k(T)^t^^^,+q(0)V',5,k'(0)) (6) 



Due to the SU (2) spin rotation symmetry, the suscepti- 
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bility function has the following form: 



^Xls?CT0Q • (?75 + ^XO^st^PaS^S (7) 



where xi ^iid xo correspond to the correlation func- 
tions of the triplet fields (such as spin) and the sin- 
glet fields (such as charge density), respectively. All the 
physical susceptibilities are determined by some compo- 
nents of Xo,ist'- For example, the total spin susceptibil- 
ity is given by xs = 5 X^s p Xiff- At the one-loop level, 
we have xosK^l;*^) = XIsKQj*^); which we denote by 
Xst'lq,^^)- For a given (q,«rj), Xol? and xu! are 4 x 4 
matrices, and the RPA susceptibility is obtained from the 
matrix equation 



xJ^(^f(q,»f^) = x(q,«f^)(ii-7o(i)x(q,*f^)) ' (8) 



with 



71 = Uhxi, 70 



U 



u 



\ 



(9) 



in the basis (si) = (11,21, 12,22). The one-loop suscep- 
tibility XsKq,*^^) is given by 



xlK^:,^^) = 



p+ q 



(27r)2 

'T-_F(-Ei/,k+q 



V -E,y,k+q — £^;y',k 
n_F(^^l/',k)) 



(10) 



with al -^ defined by Eq. 

In Fig. [2] (a), the one loop spin susceptibiHty xs(q, ^ = 
0) versus momentum q for /i = 2.0 is shown as the solid 
curve. The dashed curve shows the maximal eigenvalue 
of the one-loop susceptibility matrix x^? along the same 
contour. From this, we see that there is a critical value 
C/c — 3, at which the S* = 1 generalized RPA suscepti- 
bility diverges at an incommensurate wave vector near 
q ~ (7r/2,7r/2). This divergence occurs in the spin-one 
part of the particle-hole channel and reflects a superposi- 
tion of particle-hole spin-one fluctuations involving both 
orbitals. The RPA spin susceptibility for U — 2.8 is 
also shown in Fig. [2] (b), which, as expected, shows the 
strongest enhancement near q= (7r/2,7r/2). 

Superconductivity - Within an RPA approximation, 
the singlet and triplet pairing vertices are given bv|25| 
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To^Jik,k',^n) = -^(C/o-3C/i)*',(k-k',zf}) 



ps 



T,^J{k,k',^n) = --iUo + U,)ll{k~k',ti}) (11) 

= [ho ^ 



Here Uo^ 



tq 



7oXr^7o];: and U, 
7ixf^^7i]ps describe the effective interaction 
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FIG. 2: (a) The static one loop spin susceptibility (solid line) 
and the largest eigenvalue of one loop susceptibility matrix 
x'st along the (0, 0) (tt, 0) — > (tt, tt) (0, 0) contour in the 
Brillouin zone. The horizontal dashed line shows the value of 
1/Uc = 1/3 which indicates the critical Uc = 3. (b) The RPA 
spin susceptibility for U = 2.8. 



should be noticed that the order of orbital indices is dif- 
ferent for Uq^i and Fq,!. 

Just as for the traditional phonon case, retardation is 
important and what enters in characterizing the strength 
of the pairing interaction is 



Re[Fo(i)^«(fc,fe',c. = 0)] 



(12) 

in which a Wick rotation ifl — > cu+iS has been performed 
on Fo,i(k, k', ifi). The interaction induces scattering of 
two Cooper pairs around the Fermi surfaces. For later 
convenience, we define Ci,i = 1, ...,4 as the four pieces 
of Fermi pockets ai, a2, /32, and 7io.k the annihilation 
operator of the electron around the i^^ Fermi surface 
pocket with wavevector k G C^. Thus jia-]i is equal to 
7i^i,(Tk defined in Eq. ^ with i^i = -|-1(— 1) when Ci is 
an electron (hole) pocket. The Cooper pair defined here 
can be cither a singlet or a triplet. Here and below we 
omit the spin indices since the spin state is determined by 
the parity of the gap when k goes to — k. The scattering 
of a Cooper pair from (k, — k) on the ith Fermi surface to 
(k', — k') on the jth Fermi surface is determined by the 
projection of the interaction vertex ro,i^^(k, k') to the 
energy eigenstates: 

ro,i,^.(k,k')= J2 <,-k<kro,i^,'(fc,fe')<.k<.-lci3) 

s,t,p,q 

with k G C„ k' G Cj. 

For a pairing configuration mediated by A(k) = 
5(k)7i^_k7i,k, k G Ci, a dimensionless coupling strength 
functional is defined as[26j| 

in which i;(k) = |Vk£'i/(i)(k)| for k G is the fermi 
velocity, and <j>^ is a loop integral around the Ci 



[271 ■ ,.jps 

mediated by orbital and spin fluctuation respectively. It fermi surfaces. T\^'{k,k') = ro(i) (k, k') when ^(k) has 



4 



even (odd) parity, respectively. For a given Fq,! (k, k'), 
the optimum pairing configuration and corresponding A 
can be determined by solving an eigenvalue problem 



E 



dk' 



-^Jc, (2^)2«(k') 



-Fif'(k,k')5(k')-A.g(k), (15) 



which is obtained from the stationary condition 
SX[g{k)]/Sgik) = 0. 
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FIG. 3: (a) The effective interaction ro++(k, k') describing 
the scattering of singlet Cooper pairs on and between electron 
pockets. One wavevector k is fixed and the color shows the 
value of ro++(k, k') as a function of k'. (b) The same plot 
as (a) for the triplet channel ri(k, k'). (c) and (d) shows 
the optimum singlet and triplet pairing configurations. The 
arrow in each figure indicates a typical inter Fermi surface 
scattering process. 

The interaction Fo,iy(k, k') contains various intra 
Fermi surface and inter Fermi surface scattering pro- 
cesses. However, for /i = 2.0 we find that scattering 
on and between the f3\ and fii electron pockets is dom- 
inant. From the definition we see that Fo,ij^ only 
depends on the band label Vj. Thus the scattering on 
and between the two electron pockets are determined by 
ro,i++(k, k') Eq. p3|) . since Vi = Vj = (which we 
have shortened to -I-). The distribution of the effective 
singlet ro++(fc,A:') and triplet Fi _|__|_(fc. A;') interactions 
for V = 2.8 are shown in Fig. [3] (a) and (b) respectively 
for a fixed k G /32- From the figure it can be seen that the 
interaction for singlet pairing is repulsive, which favors a 
pairing configuration with a node. The interaction in the 
triplet channel has a smaller magnitude than the singlet, 
but is still sizable and can support a p-wave triplet state. 
The relative sign of the order parameter on the two fermi 
pockets is determined by the inter Fermi surface scatter- 
ing. For example, Fo_|__|_(k, k') is repulsive for the choice 
of k and k' shown by the arrow in Fig. [3] (a) and (c), so 
that the pairing amplitude has opposite signs at k and 
k'. For a similar reason the pairing amplitude in Fig. [3] 
(b) and (d) has the same sign at k and k". 



By solving Eq. (|15p . we obtain the optimum singlet 
and triplet pairing configurations shown in Fig. [3] (c) 
and (d) respectively. We find that Aq = 0.46 for singlet 
pairing and Ai — 0.20 for triplet pairing. To see the 
contributions of intra Fermi surface and inter Fermi sur- 
face scattering processes, we calculate these two terms 
separately by defining A^j as the term in Eq. p4)) that 
involves the scattering from the i to the j Fermi surface. 
The total A is given by A = X)i j ^ij- For i,j — /3i,/?2, 



Ao, 



are 2x2 matrices, which for t/ = = 2.8 are 



An 



0.15 
0.083 



0.083 
0.15 



Ai 



-0.026 
0.10 



0.10 
0.025 



(16) 



Here one see that the inter Fermi surface scattering makes 
an important contribution to the pairing strength A. It 
should be noticed that the singlet pairing configuration 
has the same diagonal term Xa for i = (3i and (32 while 
the p-wave diagonal terms Xa can have different values. 
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FIG. 4: The singlet Ao and triplet Ai pairing strength as a 
function of U. 

We have also studied the dependence of the pairing 
strength upon the interaction U. As shown in Fig. 21 
the strength of both the singlet d-wave and triplet p-wave 
pairing channels is increased as U increases. Within the 
RPA approximation, the coupling strength formally di- 
verge as U approaches the critical point Uc ~ 3.0 associ- 
ated with the onset of order in the spin one particle-hole 
channel. This is similar to the behavior found for the 
dx^-y2 coupling in the two-dimensional Hubbard model 
when an RPA approximation is used to treat the pair- 
ing due to the exchange coupling of spin-flucturations. 
Just as for that case, one needs to go beyond the RPA 
to determine the actual behavior of the model. 

In order to gain further insight into the nature of 
the pairing, it is useful to determine the real space 
pairing structure which corresopnds to A(k) shown in 
Fig. [3] (c) and (d). This can be obtained from A^^ — 

Sfc 5(^)7+1 (^)t+x(~^) transforming the band oper- 
tor to orbital operators, Eq ([3]), and Fourier transforming 
to the lattice coordinates. As discussed earlier, the states 
associated with the two electron pockets /?i and (i2 are 
formed primarily from xz and yz orbitals, respectively, 
as shown in Fig. [T] (b). With this in mind, we find that 
the singlet pairing in Fig 5(c) corresponds to a superpo- 
sition of singlets formed from electrons in near-neighbor 
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dxz orbitals mimus a similar super position involving the 



dyz orbitals. 



N 

r,S—x,y 

/,t 



(17) 



Under a 90 degree rotation ipj^^ and ipj^^ 

so that Ajj changes sign, corresponding to a d-wave gap. 

For the p-wave triplet shown in Fig 5(d), we find that 

r,5—x.y 

Note that this triplet gap is associated with a near neigh- 
bor intra-orbital pairing rather than the onsite inter- 



orbital triplet pairing proposed in Ref. [21 . 

Conclusion - We have studied the pairing interaction 
associated with the exchange of particle-hole fluctuation 
for a two-orbital dxz-dyz Hubbard model. By adjusting [9 
the tight binding parameters, one can obtain Fermi sur- 
face with hole and electron pockets which are similar to 
those found in bandstructure calculations for LaOFeAs. 
For a filling of two electrons per site, the signs of the hole 
and electron pockets are similar and the RPA spin sus- 
ceptibility becomes singular as the on site intra and inter 
Coulomb interaction U increase. This SDW singularity 
occurs at a wave vector q = (tt, 0) and (0, tt) associated 
with the nesting of the hole and electron pockets. When 
this model is doped, the hole pockets shrink and the elec- 
tron pockets become dominant. In this case, the SDW 
q — (tt, 0) singularity in the spin susceptibility is sup- 
pressed and there is a strong response in the spin one 
particle-hole channel near q k, (7r/2, 7r/2). The pairing 
interaction associated with the exchange of these fiuctu- 
ation leads to an attractive interaction for both singlet d- 
wave and triplet p-wave pairing, which compete closely. 
The singlet d-wave pairing strength grows faster than 
the triplet p-wave pairing strength as the interactions 
are increased and a magnetic instability is approached. 
However, more refined numerical calculations beyond the 
RPA approximation are needed to uniquely select among 
the two competing pairing states. 
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Note added - After completing this work, we learned 
that a similar work has been done by Z.-J. Yao, J.-X. Li 



and Z. D. Wang, which is posted in arXiv:0804.41 16[27|. 
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